Mesoscopic effective material parameters for 
thin layers modeled as single and double grids 
of interacting loaded wires 

Pekka Ikonen 1 , Elena Saenz 2 ' 1 , Ramon Gonzalo 2 , 
Constantin Simovski 1 , and Sergei Tretyakov 1 

1 Radio Laboratory/SMARAD Centre of Excellence 
TKK Helsinki University of Technology 
P.O. Box 3000, FI-02015 TKK, Finland 

2 Department of Electrical and Electronic Engineering 

Public University of Navarra 
Campus Arrosadia E-31006, Pamplona-Navarra, Spain 

Address for correspondence: 

Sergei Tretyakov 
Radio Laboratory, TKK 
P.O. Box 3000, FI-02015 TKK, Finland 

Fax: +358-9-451-2152 

E-mail: sergei.tretyakov@tkk.fi 

Abstract 

As an example of thin composite layers we consider single and double grids 
of periodically arranged interacting wires loaded with a certain distributed 
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reactive impedance. Currents induced to the wires by a normally incident 
plane wave are rigorously calculated and the corresponding dipole moment 
densities are determined. Using this data and the averaged fields we assign 
mesoscopic material parameters for the proposed grid structures. These pa- 
rameters depend on the number of grids, and measure the averaged induced 
polarizations. It is demonstrated that properly loaded double grids possess 
polarization response that over some frequency range can be described by 
assigning negative values for the mesoscopic parameters. Discussion is con- 
ducted on the physical meaningfulness to assign such material parameters 
for thin composite slabs. The results predicted by the proposed method for 
the double-grid structures are compared with the results obtained using the 
commonly adopted S'-parameter retrieval procedure. 
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1 Introduction 



The present work is motivated by recent articles in which the authors assign 
very exotic material parameters e, fi (or effective refractive index n) for thin 
composite slabs consisting only of a few grids of interacting complex-shaped 
inclusions, e.g. [1-4]. In all of the above referred works the effective material 
parameters have been extracted using the ^-parameter retrieval procedure 
whose underlying basics are described e.g. in [5-7]. In this procedure one 
replaces the actual structure by a slab of homogeneous material having a 
certain effective thickness and searches for effective material parameters that 
will reproduce the scattering results obtained with a normal plane-wave in- 
cidence. The method is often known to lead to antiresonant behavior for 
one of the extracted material parameters, which further leads to physically 
wrong sign for the imaginary part of the corresponding material parame- 
ter [8]. According to the classical electrodynamics [9], a physically meaning- 
ful material parameter is a linear response that must be independent of the 
electromagnetic field phasors excited in the medium, and especially in case of 
composite media, independent of the geometrical construction constituting 
the medium. Moreover, the causality requirement (which is directly related 
to the fulfillment of the Kramers- Kroning relations) and the passivity re- 
quirement (which relates to the allowed sign of the imaginary part) must be 
simultaneously fulfilled [9-12]. It is therefore clear that an effective parame- 
ter showing antiresonant behavior does not bear the meaning of a physically 
sound material parameter describing a real homogeneous materiatj, it only 
allows to reproduce the scattering scenario with a specific excitation [13]. 
When adopting the aforementioned classical definition for physically sound 
material parameters it is clear that, according to this definition, the param- 
eters presented e.g. in [1,2,4] are nonphysical. Also, in the view of classical 
homogenization approaches (e.g. [14]) it is clear that a structure consisting 



1 Please note that the function typically fitting into the observed antiresonant frequency 
dependence [8] resembles the Lorentz model with the sign in front of the fractional term 
changed (we assume e Jut time dependence): 

F( U ) = 1 A ^ 



uJq - uj 2 + jtoT 



The change of the sign assures that the poles of the function lie in the upper half-omega 
plane (thus the function can be said to be causal), however, the sign of the real part of 
frequency derivative is wrong to fulfill the passivity requirement. 
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only of a few (e.g. two) interacting grids of inclusions does not correspond 
to any slab of an effectively homogeneous material. In an effectively homo- 
geneous material the unit cell should contain a large number of molecules 
or inclusions, and yet still remain small compared to the wavelength. This 
partly assures e.g. that the change of the material sample thickness does not 
alter the polarization response of a small material volume. Evidently, the ma- 
terial parameters assigned for thin composite slabs are mesoscopic, thus, they 
depend on the number of monolayers (grids) constituting the sample. Most 
metamaterials and thin composite layers operate in the mesoscopic regime, 
intermediate between the homogeneous medium regime, and the regime of a 
collection of weakly interacting individual scatterers. 

Typically the authors using the ^-parameter retrieval procedure explain 
that the finite spatial periodicity of the sample causes the often observed 
antiresonant behavior. It is also well known that the method introduces an 
ambiguity in the phase of the wave reflected from finite thickness structures. 
To avoid this ambiguity, attempts have been made to assign effective mate- 
rial parameters for finite thickness complex slabs through refraction, rather 
than through reflection. For example, in [15] the authors a present clear dis- 
cussion on the challenges to assign physically sound effective refractive index 
for wire medium slabs having a finite thickness. A small artificial interface 
shift is introduced for the actual structure and the authors use an iterative 
minimization technique to find the effective refractive index of a homoge- 
neous slab that transmits the same amount of light as the actual complex 
structure (essentially the same technique was used also e.g. in [16]). Authors 
of [15] also define the refractive index through reflection and transmission, 
but they search for the material parameters that allow one to reproduce the 
actual scattering results with all incidence angles (note that the technique in- 
troduced e.g. in [5] uses only reflection and transmission data for a normally 
incident plane wave). 

Motivated by the recent contradictory dispersion results obtained for sin- 
gle and double grids of interacting inclusions [1,2,4] we turn our attention 
to single and double grids of loaded wires. Grids of loaded wires are con- 
sidered since the theory of interacting grids of such wires is well established 
and accurate [17,18], and certain specific loading scenarios allow us to sim- 
ulate the dispersion properties of lattices of small dipole scatterers (more 
discussion on this will be presented later). Also, despite the fact that the 
local field approach used in this paper completely takes into account the 
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near-field interactions (all the higher-order Floquet harmonics) between the 
wires, is remains physically illustrative and computationally efficient. Our 
goal is to assign mesoscopic material parameters for the proposed grids of 
wires directly by solving the induced dipole moment densities and averaged 
fields. Recently we used the same approach for assigning effective material 
parameters for grids of interacting inclusions modeled as small electric dipole 
scatterers [19]. In this work our goal is, in addition to studying different 
structures as considered in [19], present a more extensive discussion on the 
physics behind the obtained polarization responses, and more systematically 
compare (with double grid structures) the proposed results to the results 
obtained using the ^-parameter retrieval procedure. The details of deriva- 
tion are introduced, and the obtained results and their physical meaning are 
thoroughly discussed. 



2 Formulation 

2.1 Outline of the strategy 

The known constitutive relations in a certain macroscopic medium (or in a 
large sample of a certain macroscopic material) read [20] 

D = e E + P, (1) 

B = /i H + M, (2) 

where D and B are the averaged electric and magnetic flux density, E and 
H are the averaged electric and magnetic field, P is the averaged polariza- 
tion (electric dipole moment density), and M is the averaged magnetization 
(magnetic dipole moment density). When introducing a physically small 
averaging volume V, the average polarization and magnetization can be cal- 
culated as [20] 

p = F = ^/ i(r '>" v ' (3) 

M = V = vI T ' x i(r,) dV ' (4) 

where i is the microscopic current density confined in the averaging volume, 
and the integration is done over that volume. 



5 



In the following analysis our strategy is to first solve the current induced 
to the wires by the incident field, taking into account the full interaction 
between all the wires. After this, we define an artificial unit volume and 
calculate the macroscopic polarization and magnetization using eqs. ([3]) and 
(H]). Further, the total averaged fields in the unit volume are defined, and 
eqs. ([T]) and (j2J) are used to express the mesoscopic material parameters in 
terms of macroscopic polarization (magnetization) and averaged fields. 



2.2 Single grid of loaded wires 

Consider a periodical grid of identically loaded, infinitely long wires depicted 
schematically in Fig. [TJ The radius of wires r$ is assumed to be small enough 
so that the wires can be represented as infinitely long current lines located at 
the position of the wires axis. Throughout the rest of the paper we consider 




Figure 1: Grid of thin wires loaded with bulk impedances (Z = corresponds 
to a grid on unloaded wires). The wires are infinitely long and the number 
of wires in the y-direction is infinite. The insertion period for the loads I is 
much smaller than the wavelength. An example unit cell is shown shadowed. 
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normal plane-wave incidence and express the external electric field as 



E cxt = u z E cxt e~ jkx . (5) 

Let us first assume the wires to be unloaded and pick up a reference wire 
from the grid. The relation between the amplitude of the local electric field 
E l ° c acting on the surface of the reference wire and current / flowing on the 
wire surface reads 

E l ; c = ao 1 /, (6) 

where ao is the susceptibility of the wire which can be determined from the 
boundary condition on the wire surface (see e.g. [23]): 



a 



(7) 



(2) 

where rj is the wave impedance in the background medium, and Hq rep- 
resents the Hankel function of the second kind and zero order. If the wires 
are loaded with a certain uniformly distributed impedance per unit length 
Z, the inverse of the wire susceptibility becomes [18,22]: 

a' 1 = a 1 + Z. (8) 

In the following we denote the wire susceptibility generally as a (qzq is a 
special case for the susceptibility when Z = 0). The amplitude of the local 
electric field acting on the reference wire can also be expressed in the following 
manner: 

where E ey± is the amplitude of the external field at the surface of the reference 
wire and (3(0) is so called self interaction coefficient given by eq. (11) in [17] 
(see also [21]): 



m = -f 



1 ,j ( \ kd \ j ^ ( 1 d 



kd 2 7T V 4tt 7 d^ Q \ ^/(27in/d) 2 -k 2 27r|n 

(10) 

where d is the separation of the wires in the grid, and 7 ~ 0.5772 is the 
Euler constant. Physically (3(0) tells how strongly the wires located in the 
same grid as the reference wire contribute to the amplitude of the local field. 
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From eqs. and §B) we can solve the amplitude of current induced to the 
reference wire as 

^cxt 

Next, let us turn our attention to the definition of the material parameters. 
Let the unit cell volume (xxyx z) shown shadowed in Fig. \l}beV = sxdxl. 
When the integration in eq. ([3]) is conducted, the average polarization reads: 



P = u 2 — = — , (12) 
jujsa jus 

where we have defined the unit cell averaged surface current density J = 
u z I/d. We proceed on defining the total averaged electric^] field in the grid 
plane. In the far zone, the electric field scattered by the grid (in the normal 
direction) is a plane-wave field, and close to the grids the amplitude of this 
field can be expressed as [22] 

E sc = (13) 

The far field reflection and transmission coefficients (referred to the grid 
plane) can be formulated using the ratio between the scattered and external 
electric fields [22] 

rpsc rpsc 

R = —r, T=l + —. (14) 

By definition [22, p. 79], the total averaged electric field in the grid plane is 
the sum of the external field and the plane wave field created by the grid: 

E = E cxt - |j. (15) 

To calculate the averaged field we need to average E over the unit cell thick- 
ness s (since we consider plane waves the only variation in the fields occurs 
in x-direction). The final result reads: 

g = sin(fa/2) F , _ (16) 



ks/2 jks/2 2 



2 In the case of a single grid of wires we do not need to consider the averaged magnetic 
field, as due to the thin wire assumption there is no circulating current, and the amplitude 
of the induced magnetic dipole moment density is zero. 



8 



When the unit cell thickness s tends to zero the limit obviously reduces 
to eq. (TT51) . By inserting eqs. (lT2"j) and (1TE1) into and dropping out the 
vector notation we get the following expression for the mesoscopic effective 
permittivity: 

kJ , , 

e = l + ^ — . 17 

juje [2sm(ks/2)E cxt + jr}J(l - e~i ks / 2 )] 

When we use relation ffTTT) to express J as a function of E ext another expres- 
sion is obtained: 

(18) 



2d/r j [a~ 1 - (3(0)} sin(fcs/2) + e~i ks / 2 - 1 ' 

One should bear in mind that the current induced to the grid varies as 
a function of the incidence angle, thus, different polarization response is ex- 
pected to occur when the incidence angle varies [21,22]. Importantly we also 
observe that for the oblique incidence also higher order Floquet harmonics 
contribute to the scattered electric field calculated over the unit volume thick- 
ness. Neglecting the effect of these harmonics would leave us with non-local 
mesoscopic permittivity [12]. 

The fact that the above deduced permittivity is not independent of the 
field phasors excited in the slab (incidence angle), or the geometrical con- 
struction (thickness) constituting the slab clearly indicates that the permit- 
tivity function does not fulfill the classical criteria set for physically sound 
material parameters of a real homogeneous material [9]. Nevertheless, when 
accepting this fact and bearing in mind the mesoscopic nature of the deduced 
parameter, the permittivity expression is still useful in studying practically 
engineered slabs (although they are not sufficient for enough complete de- 
scription of the composite): correctly deduced mesoscopic material param- 
eters contain physically sound information about the actual (microscopic) 
polarization response of the structure. 



2.3 Double grid of loaded wires 

In this subsection we consider two grids of wires as shown in Fig. [2j In this 
case at some frequency (and with proper impedance loading) currents in the 
grids will be out of phase giving rise to a rather strong circulating current and 
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induced magnetic dipole moment density (this will be demonstrated using a 
specific example later). Following the reasoning in the previous subsection we 
first write down the system of equations for the unknown induced currents: 

a- l h = E^e ]kh + (3(0)1! + (3(2h)I 2 , 

a~ l I 2 = E cxt e~' kh + (3{2h)h + (3(0)I 2 , (19) 

where (3(2h) is so-called mutual interaction coefficient (eq. (12) in [17]) taking 
physically into account the mutual interaction between the grids: 

^ = - y y t—r^. (20) 

1 K ' 2kd ^ k (n) 

n=-oo ^ 

Above k^ is the ^-component of the n-th Floquet mode wave vector [17]: 



fcjW = -jJ(ky + 27rn/dy + k 2 z - k 2 , Re{v^)} > 0. (21) 



s 




Figure 2: Two grids of wires loaded with bulk impedances (Z = corresponds 
to grids of unloaded wires). The wires are infinitely long and the number of 
wires in the y-direction is infinite. The insertion period for the loads I is 
much smaller than the wavelength. An example unit cell is shown shadowed. 
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Currents in the grids can be found by solving the system of equations (I19p . 
The average polarization and magnetization take the following forms after 
applying eqs. (j3J) and (TjJ (we use the notation introduced in the previous 
subchapter for the unit cell dimensions): 



Ji + J 2 

jus 



(22) 



M = /i -(J 2 

s 



Ji). 



(23) 



Both grids will create a scattered electric field whose amplitude is given 
by eq. ffl3|) . When the electric field is averaged over the unit cell thickness 
the result reads: 



E 



E ext sin(A;s/2) + + J 2 )(l - cos(kh)e 



-jks/2\ 



ks/2' 



(24) 



The expression for the averaged magnetic field reads: 



V 



E cxt sin(fcs/2) - ^(Ji - J2) sm(kh)e 



-jks/2 



_ ks/2' 



(25) 



Further, the mesoscopic effective material parameters are determined by ap- 
plying eqs. ([1]) and ([2]). The final expressions read (in the double-grid case 
the expressions become quite cumbersome when the current densities are sub- 
stituted as functions of the external field. For this reason only expressions 
containing the current densities are presented): 



1 + 



k(Ji + J 2 



jLue {2E cxt sin(ks/2) + jt]{Ji + J 2 )[l - cos(kh)e-i ks / 2 }} 

khr)(J 2 - J\) 



/i = 1 



2E ext sin(fcs/2) - //(Ji - J 2 ) sm(kh)e-i ks / 2 



(26) 
(27) 



3 Example results 

Let us choose the following structural parameters for the wire grids (see 
Figs. [T]and [2]): r = 0.1 mm, d = 200r . The insertion period for the loads 
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I = 50r . In the double-grid case the separation between the grids 2h = 40r . 
These parameters are kept the same throughout the rest of the paper. 

We consider three different loading scenarios: 

1. the wires in the grid(s) are unloaded (Z = 0), 

2. the wires are loaded with bulk capacitors [Z = l/(ju>Cl)], 

3. the wires are loaded with a parallel connection of bulk capacitors and 
inductors [l/Z = 1(jluC + l/{jujL))- 1 ]. 

Please note that we do not separately consider the loading of the wires with 
bulk inductors or with a series connection of bulk capacitors and inductors. 
These two situations effectively resemble the scenarios described in items 1) 
and 2), as the effect of inductive loading in series can be attributed to the 
change of the effective wire radius [18]. 

3.1 Single grid of wires 
3.1.1 Unloaded wires 

Let the wires in the grid be first unloaded [Z = 0). The reflection and 
transmission coefficients calculated using eq. (THI) are depicted as a function 
of the normalized wire separation in Fig. [3J We qualitatively repeat the well 
known observation (e.g. [21]) according to which at low frequencies the grid 
is strongly reflective (the wave does not probe the fine details of the grid). 
Note that the normalized frequency d/X — 1 corresponds to a certain kind 
of resonance [21,22]: at this frequency the contribution from the other wires 
to the local electric field acting on the surface of the reference wire adds up 
in phase. However, the field amplitude at the wire surface must be limited 
which is only possible when the current induced to the grid is zero, thus, the 
grid is ideally transparent at this frequency. 

For the calculation of the mesoscopic permittivity we need to choose the 
unit cell thickness s. One rather intuitive choice is s = 2r . In addition to 
this, we partly follow the analysis of [15, 16] where the the authors use an 
artificial interface shift d/2 on both sides of stacks of wire grids when repre- 
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senting the actual structures as slabs made of homogeneous materials Thus, 
in the second scenario s = d (when representing the grid as a homogeneous 
slab we imagine that there is an artificial interface shift d/2 on both sides of 
the grid). 

The permittivities for the both aforementioned situations, calculated using 
eq. (fl8|) . are shown in Fig. HI We notice that in both cases the permittivity 
is purely real as it should be since the wires (and the background medium) 
are lossless. In the case s = 2r$ the frequency derivative of permittivity 
is positive over the entire frequency range, thus, the permittivity function 
strictly satisfies the passivity requirement. Interestingly, however, this is not 
the case when we introduce the artificial interface shift: the permittivity 
has a maximum roughly at d/X = 0.72, and the slope becomes negative 
as the frequency approaches the point d/X — 1. At these frequencies the 
permittivity function will lead to a physically wrong sign for the imaginary 
part when inserted into the Kramers-Kroning relations. However, also the 
wire spacing is comparable to the wavelength at these frequencies and the 
validity of the effective medium description is weak. At frequencies below 

3 In [15,16] the lattice constant is the same in both transversal (with respect to the wire 
axis) directions. 
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Figure 3: Reflection and transmission coefficient for a single grid of unloaded 
wires as a function of the grid's spatial frequency. 
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Figure 4: Mesoscopic effective permittivity of a single grid of unloaded wires 
as a function of the grid's spatial frequency. 

d/X « 0.72 we can, however, at least qualitatively make some interesting 
observations from the results depicted in Fig. HI 

Typically in the quasi-static regime the effective permittivity of infinite 
lattices of unloaded wires (when excited with normally incident planes waves) 
is described using the lossless Drude model (e.g. [17,24]): 



where c^o is the angular "plasma" frequency which is determined by the 
structural parameters of the grid. If we assume that instead of a single grid 
of wires we would investigate the properties of an infinite wire lattice with 
rectangular periodicity (the lattice constant in both transversal directions 
would be d), the normalized wire medium plasma frequency, estimated us- 
ing equations available in [17,24], would be d/\ Q = 0.20. When s = d we 
estimate from Fig. H]that d/Xo ~ 0.20. Actually, a similar observation has 
been made already in [16]: the authors observed that in case of an interface 
shift corresponding roughly to d/2 the effective permittivity of an equivalent 
homogeneous slab simulating the single grid closely approximates the per- 
mittivity of an infinite lattice. This is intuitively rather clear after simple 
symmetry considerations [24]. When the unit cell thickness s = 2ro the be- 
havior of the permittivity function is dramatically different compared to the 



e = 




(28) 
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Figure 5: Reflection and transmission data calculated for the homogenized 
slabs simulating a grid of unloaded wires. 

case s = d. 

Next, we represent the actual wire grid as a slab of homogeneous material 
having thickness (in x-direction, Fig. [T|) 2r or d. The corresponding effec- 
tive permittivities characterizing these homogeneous slabs are those shown in 
Fig. HI We use standard equations for calculating the reflection and transmis- 
sion coefficients through such homogeneous slabs, see e.g. the transmission 
matrix formalism presented in [22, Ch. 2.4]. The results are depicted in 
Fig. El 

In the first case the effective slab thickness equals to the diameter of the 
wires, and we observe a very close agreement between the exact result de- 
picted in Fig. [31 and the result calculated for the homogeneous slab. For 
this specific example we have, thus, verified that the effective permittivity 
depicted in Fig. @] correctly describes the actual (microscopic) polarization 
response of the grid, and also allows to reproduce the reflection and trans- 
mission results with a normal plane- wave incidence. 

When the slab thickness equals the wire spacing the agreement between 
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Figure 6: Reflection and transmission coefficient for a single grid of wires 
loaded with bulk capacitors as a function of the grid's spatial frequency. 

the exact result and the result depicted in Fig. is good roughly up to 
d/X — 0.1. At higher frequencies the discrepancy is visible (please remember 
also the discussion concerning the meaningfulness of the permittivity function 
for this case at frequencies above d/X = 0.72). It is evident from the results 
that the effective permittivity for this situation (see Fig. @| does not allow to 
reproduce the actual reflection and transmission results with normal plane 
wave incidence. Obviously the difficulties arise from the correct definition 
of the equivalent homogeneous slab thickness. However, permittivity defined 
when s = d gives a good estimate for the macroscopic polarization response of 
an infinite wire lattice (at least up to frequencies d/X = 0.72, in our example 
case) . 

3.1.2 Capacitively loaded wires 

Let us next load the wires by bulk capacitors whose value is chosen to be 
C = 1 pF. It is known [18] that when the load capacitance value is infinitely 
small, the dispersion characteristics of a lattice of such wires resemble the 
characteristics of a lattice of point dipole scatterers. When a realistic load 
capacitance value and the insertion period are correctly chosen, this equiv- 
alence holds still, to some extend, at low frequencies [18]. Thus, one can 
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consider a properly capacitively loaded wire grid as a grid of short wire 
(dipole) scatterers. Non-physical dispersion results (defined through reflec- 
tion and transmission) for grids of such inclusions have been presented earlier 
e.g. in [2,4]. 

The reflection and transmission coefficients calculated using eq. (fT4|) are 
shown in Fig. [16j We see that there appears a stop-band not observed in the 
previous case (grid of unloaded wires). The behavior seen in Fig. (THlcan be 
explained (qualitatively) by regarding the grid as a bulk load, consisting of a 
series connection of inductance and capacitance, in a waveguide representing 
free space wave propagation [22]. At low frequencies the impedance of the 
grid is strongly capacitive and the load represents an open circuit in the 
waveguide, thus, waves can propagate through the grid. At some frequency 
the wires are in self-resonance, and in this case the waveguide is effectively 
short-circuited leading to a stop-band [18,22]. 

The mesoscopic effective permittivity calculated using eq. (fI5|) for this 
case is shown in Fig. The permittivity is purely real, and we observe 
that it obeys a Lorentz-type behavior, which is consistent with the effective 
permittivity deduced for a semi-infinite capacitively loaded wire medium [18]: 

c=l + -5-^, (29) 



UJq — UJ 2 

where A is an amplitude factor depending on the load capacitance value and 
the structural parameters of the grid (also uq in this case depends on both 
the structural parameters and the load capacitance value). This behavior is 
physically due to the fact that in the self-resonance the distributed impedance 
value on the wire surfaces is low and allows strong current to be induced to 
the wire. Strong induced current directly means strong polarization response, 
therefore the very high (resonant) value of the permittivity. Moreover, the 
choice of the unit cell thickness s again dictates the strength of the polariza- 
tion response. Please note that in reality, due to finite conductivity of the 
wires and parasitic losses in the lumped elements, the effective permittivity 
would show a significant imaginary part in the vicinity of the resonance as 
is imposed by the Kramers-Kroning relations. Similarly as with the grid of 
unloaded wires we observe that the permittivity, when defined using s = 2ro, 
is a monotonically growing function outside the resonant band. A closer 
investigation of the permittivity function, defined when s — d, reveals that 
the slope of the function is negative between the frequencies d/X = 0.72 ... 
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(a) Permittivity over a wide frequency range. 



(b) The same result as in (a), but plotted over a 
more narrow frequency range. 



Figure 7: Mesoscopic effective permittivity of a single grid of wires loaded 
with bulk capacitors as a function of the grid's spatial frequency. 



d/X = 1. The behavior at these frequencies is therefore similarly nonphys- 
ical as with the grid of unloaded wires (however, also the effective medium 
description is weak). 

Authors of [24] considered an infinite lattice of capacitively loaded wires 
and deduced the following quasi-static limit for the effective permittivity 
in case of relatively high capacitive impedance (the notation is clear from 
Fig. CD : 

ni 

e s = 1 + —j. (30) 
e sd 

2r and s 



When the unit cell thickness is s = 2r and s = d, eq. fl30|) predicts e s = 142 
and e s = 2.4, and the the corresponding values read from Fig. 7(a) are exactly 
the same. For an infinite lattice of capacitively loaded wires with d as the 
lattice constant in both transversal directions and the structural parameters 
considered here, the normalized plasma frequency is roughly d/Xo = 0.168 
[18]. From Fig. Owe see that the corresponding frequencies are d/X^ = 0.180 
(s = 2r ) and d/X = 0.163 (s = d). Thus, again the introduction of the 
artificial interface shift allows us to rather well approximate the refractive 
properties of an infinite lattice (at least in terms of the effective plasma 
frequency and static effective permittivity). 
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(a) Unit cell thickness s = 2tq. 
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(b) Unit cell thickness s = d. 



Figure 8: Reflection and transmission data calculated for the homogenized 
slabs simulating a grid of capacitively loaded wires. 

The reflection and transmission coefficients calculated for the homoge- 
neous slabs are shown in Fig. [HI We observe that when the unit cell thickness 
equals the wire diameter the mesoscopic effective permittivity, depicted in 
Fig. [7J represents the grid correctly both from the polarization and reflec- 
tion/transmission points of view. If an artificial increase in the slab thickness 
is introduced the reflective/transmissive properties dramatically differ from 
the exact results close to the self-resonant frequency (as with the grid of un- 
loaded wires the results agree well up to d/ X — 0.1). Near that self-resonance 
the absolute value of the effective refractive index is very large, thus, the 
electrical slab length is very large. In a slab whose physical thickness no- 
ticeably differs from zero such condition can be predicted to cause thickness 
resonances that would be seen as sharp transmission maxima. This, indeed 
seems to be the case according to Fig. [HJ 



3.1.3 Wires loaded with a parallel connection of capacitance and 
inductance 

In this case the wires are loaded with bulk impedances consisting of a parallel 
connection of capacitance and inductance. The capacitance value is kept the 
same as before, the load inductance value is L — 1 nH, and there appears a 
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parallel resonance at d/\\ = 0.34. The reflection and transmission coefficients 
for such a loaded grid are depicted in Fig. O Once again, the behavior can be 
straightforwardly explained using the idea of a waveguide loaded (in parallel) 
with a bulk load. At low frequencies the the impedance of the load is strongly 
inductive and the waveguide is effectively short circuited. When the load 
circuitry is in parallel resonance the impedance of the load tends ideally to 
infinity, and waves can propagate through the grid. With our load values and 
structural parameters a series resonance occurs just above the parallel load 
resonance. This resonance is caused by the inductive wire segments (that 
can be modeled as bulk inductors) and load capacitors, and is observed as 
the sharp reflection maximum in Fig. O 

When looking at the calculated mesoscopic effective permittivity, shown 
in Fig. [TUl we observe that the permittivity crosses unity at the parallel load 
resonance (this is seen as the crossing point of the permittivity functions in 
Fig. [TO]) . This explains also the total transmission at this frequency. Around 
the frequency point of series load resonance the behavior is qualitatively the 
same as in the case of capacitively loaded wires. Note again that the permit- 
tivity is purely real in both cases. However, again the frequency derivative 
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Figure 9: Reflection and transmission coefficient for a single grid of wires 
loaded with a parallel connection of bulk capacitors and inductors as a func- 
tion of the grid 's spatial frequency. 
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(a) Permittivity over a wide frequency range. 




s = 2r„ 
s = d 
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(b) The same result as in (a), but plotted over a 
more narrow frequency range. 



Figure 10: Mesoscopic effective permittivity for a grid of wires loaded with a 
parallel connection bulk capacitors and inductors as a function of the grid's 
spatial frequency. 




Figure 11: Reflection and transmission data calculated for the homogenized 
slabs simulating a grid of capacitively and inductively loaded wires. 
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of e, defined using s = d, is negative at frequencies higher than d/X = 0.72. 

The reflection and transmission coefficients calculated for the homoge- 
neous slabs are shown in Fig. [TTJ The behavior is qualitatively the same 
as already described with the previous examples: when the thickness of the 
homogeneous slab equals the wire diameter, the slab correctly describes the 
actual grid both from the polarization and reflection/transmission points of 
view. 

In short intermediate conclusion related to the single grids, the mesoscopic 
effective permittivity, defined directly through induced electric dipole mo- 
ment density and averaged electric field, is a physically sound function with 
all the loading scenarios over the spectral range where the effective medium 
description is valid. Moreover, the calculated examples indicate that the 
grid can accurately be represented both for polarization (equivalently, for re- 
fraction when defined through induced dipole moment density and averaged 
electric field) and reflection/transmission as a homogeneous slab having the 
proposed permittivity when the unit cell thickness equals the diameter of the 
wires. This observation holds when the diameter is very small. When an ar- 
tificial increase in the unit cell thickness corresponding to d/2 is introduced, 
the permittivity function representing the homogeneous slab does not allow 
to reproduce the reflection/transmission results for the actual grid. However, 
the function gives (up to a certain frequency) a rather good approximation 
for the effective permittivity of an infinite lattice. 

3.2 Double grid of wires 
3.2.1 Unloaded wires 

The reflection and transmission coefficients (referred to the plane of the left- 
most grid in Fig. [2]) are shown in Fig. [J2] (we have used eq. 014|) by also taking 
into account the scattered field created by the second grid, and the phase 
shift between the grids). The coefficients qualitatively resemble the results 
depicted in Fig. [3] for a single grid of unloaded wires with the exception that 
the zero point of reflection coefficient has shifted to a lower frequency due to 
the influence of the second grid. 

The mesoscopic effective material parameters of the equivalent homoge- 
neous slab are shown in Fig. [T3J In the double grid analysis we set the unit 
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cell thickness equal to s = 2h (accurately corresponding to the actual physi- 
cal thickness of the structure) and s = Ah, following the analysis of [15, 16]. 
Please note that in this case the lattice constants in the transversal directions 
(with respect to the wire axis) are different, and we set the artificial interface 
shift equal to half of the lattice constant in ^/-direction (Fig. [2]). 

We observe that all the material parameters are purely real, and the fre- 
quency derivative of every parameter is a growing function over the entire 
spectral range. The effective permittivity behaves qualitatively in a similar 
manner as introduced in connection with the single grid. However, in this 
case the double grid is very slightly diamagnetic: the static effective perme- 
abilities for the two homogeneous slabs, defined from Fig. 13(b), are 0.80 
(s = 2h) and 0.91 (s = Ah). A similar observation related to diamagnetism 
has often been observed in connection with the results obtained using the 
^-parameter retrieval method [25] . For an infinite lattice having lattice con- 
stant 2h in x-direction and d in ^/-direction, equations available in [18] predict 
that the normalized plasma frequency of such a wire medium is d/\o = 0.50. 
The normalized plasma frequencies observed in Fig. [13] correspond roughly 
to d/X « 0.60 (a = 2h) and d/X « 0.44 (s = Ah). 

Fig. dH depicts a comparison between the effective material parameters ob- 




Figure 12: Reflection and transmission coefficient for a double grid of un- 
loaded wires as a function of the grid's spatial frequency. 
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(a) Material parameters over a wide frequency (b) The same result as in (a), but plotted over a 
range. more narrow amplitude range. 



Figure 13: Mesoscopic effective permittivity and permeability of a double grid 
of unloaded wires as a function of grid's spatial frequency. 



tained using the proposed method and the S-parameter retrieval method (for 
the S'-parameter retrieval method we have used equations available in [5]). 
In this example the thickness of the equivalent homogeneous slab equals the 
physical thickness (s = 2h), and we observe that the two methods give rather 
closely agreeing results. However, the frequency derivative of permeability 
function, obtained using the S'-parameter retrieval procedure, is very slightly 
negative over the entire plotted frequency range. The static effective perme- 
ability obtained using the S-parameter retrieval procedure equals 0.87. 

Fig. [15] shows the reflection and transmission coefficients for the homo- 
geneous slabs characterized by the mesoscopic effective parameters depicted 
in Fig. O We see that with both slab thicknesses the results agree rather 
closely with the exact results depicted in Fig. [T2J Additional calculations 
(not shown) indicate that when the interface shift is increased the disagree- 
ment between the reflection/transmission results for the actual structure and 
for the homogeneous slab increases, and eventually there appears a null in 
the reflection coefficient around d/A~ 0.5. 
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Figure 14: Comparison between the mesoscopic effective material parameters 
obtained using the proposed method and the S -parameter retrieval method. 

3.2.2 Capacitively loaded wires 

The reflection and transmission coefficients for this case are depicted in 
Fig. [161 Once again the behavior resembles the behavior of a single grid 
of capacitively loaded wires with the important exception that there appears 
a narrow transmission window inside the wider stopband. This transmission 
window originates physically from the fact that the currents in the grids can- 
cel each other at this frequency (when the phase shift between the grids is 
taken into account, remember that the reflection coefficient is referred to the 
plane of the leftmost grid in Fig. [2]) . 

The calculated mesoscopic effective material parameters are depicted in 
Fig. [13 Please note that all the material parameters are purely real and 
behave in a physically sound manner (the frequency derivative is positive ev- 
erywhere outside the resonance region). Moreover, the frequency dependence 
of both parameters clearly resembles the behavior predicted by the classical 
Lorentz model. The behavior in the effective permittivity arises from the 
self-resonance of the wires, i.e. when the load circuitry is in series resonance 
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Figure 15: Reflection and transmission data calculated for the homogenized 
slabs simulating a double grid of unloaded wires. 
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Figure 16: Reflection and transmission coefficient for a double grid of wires 
loaded with bulk capacitors as a function of the grid's spatial frequency. 
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(a) Material parameters over a wide frequency (b) The same result as in (a), but plotted over a 
range. more narrow frequency range. 



Figure 17: Mesoscopic effective permittivity and permeability of a double 
grid of wires loaded with bulk capacitors as a function of the grid's spatial 
frequency. 



the current induced to the wires will be strong leading to a strong polariza- 
tion response. The resonant permeability arises from circulating out-of-phase 
currents flowing in the grids. The difference in the current amplitudes flow- 
ing in the grids is minimized at a slightly higher frequency as compared to 
the deep reflection minimum observed in Fig. [T6j Therefore the reflection 
minimum (in this particular case) does not exactly correspond to the res- 
onant point of permeability. It is also worth to observe that resonance in 
permeability is a very narrowband phenomenon compared to the resonance 
in permittivity. This is because the mutual cancelation of the grid currents 
occurs in a narrow frequency window due to very strong mutual coupling 
(very small grid separation) between the grids. 

It is also worthwhile to notice that both mesoscopic effective parameters 
are negative in a narrow frequency band just above the resonance of per- 
meability, thus, the polarization response at this frequency range can be 
described by assigning "negative effective refractive index" for the structure. 
We have to bear in mind, however, the important physical limitations out- 
lined above for the introduction of effective parameters for thin complex 
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slabs. One should also observe that the amplitudes of permittivity and per- 
meability (Fig. [T7|) are at the same (negative) level only over an extremely 
narrow frequency band. This means that the structure (with these particu- 
lar structural dimensions and load values) would be very difficult to match 
to incident radiation. When the unit cell thickness is s = 2h and s = 4h, 
eq. (15U|) predicts e s = 7.7 and e s = 4.5. The corresponding values read from 
Fig. [T7] are e s = 15.1 and e s = 8.1. Static permeabilities for both effective 
slab thicknesses equal unity. 

Comparison between the results obtained using the proposed method and 
the S-parameter retrieval method is depicted in Fig. [18j The thickness of 
the equivalent homogeneous slab equals the physical thickness (s = 2h). 
Interestingly, the effective permittivities agree rather well over the whole 
depicted frequency range (visible disagreement is observed only close to the 
permittivity resonance point). However, the effective permeabilities obtained 
using different methods behave in a visibly different manner: the predicted 
magnetic response of the double grid is much weaker according to the S'- 
parameter retrieval results. Also, it is evident that the function predicted 
by the S'-parameter retrieval method cannot be represented using the causal 
Lorentz model. 




6JX dfk 

(a) Comparing the mesoscopic effective permittiv- (b) Comparing the mesoscopic effective permeabil- 
ities, ities. 

Figure 18: Comparison between the mesoscopic effective material parameters 
obtained using the proposed method and the S -parameter retrieval method. 



28 



The reflection and transmission coefficients, calculated when the actual 
structure is represented as a slab of homogeneous material (the material 
parameters are those depicted in Fig. [T7|) . are shown in Fig. [191 We notice 
that in the vicinity of the resonant point of material parameters there appears 
several thickness resonances. At other frequencies the agreement between the 
exact result (Fig. [TBI) and the results calculated for the homogeneous slabs 
is rather good. 

We can conclude that the effective material parameters, shown in Fig. [XT], 
correctly describe what kind of polarization mechanism takes place in the mi- 
crostructural level. However, when the actual grid structure is represented 
as a homogenized slab, the proposed mesoscopic parameters do not allow 
to exactly reproduce the reflection/transmission results. Visible disagree- 
ment is observed close to the resonant points of the material parameters. 
Close to these resonant frequencies the electrical slab thickness is very large, 
and this leads to several closely located thickness resonances not observed 
with the actual grid structure. Note that when the wires are non-resonant, 
e.g. unloaded as in the previous subchapter, the agreement between the re- 
flection/transmission results for the actual structure and homogenized slabs 
is rather good over the entire spectral range. This indicates that when the 
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(a) Unit cell thickness s = 2h. 
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(b) Unit cell thickness s = Ah. 



Figure 19: Reflection and transmission data calculated for the homogenized 
slabs simulating a double grid of capacitively loaded wires. 
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slab is electrically very dense it is difficult to assign the correct equivalent slab 
thickness reproducing the scattering scenario (intuitively this is clear since 
the double-grid cannot physically correspond to a sample of any homogeneous 
material). It it important to bear in mind that the above mentioned defi- 
ciencies are avoided with the S'-parameter retrieval method since the method 
is based on the fulfillment of the actual scattering scenario. However, as is 
shown by the above example (and will further be shown below), the material 
parameter functions obtained using the S'-parameter retrieval results might 
not satisfy the causality or the passivity requirement (see also e.g. [1-4,8]). 

3.2.3 Wires loaded with a parallel connection of capacitance and 
inductance 

Finally, we proceed on studying the situation when the wires are loaded with 
bulk connections of capacitances and inductances. The reflection and trans- 
mission coefficients are shown in Fig. [20j an d they again resemble the results 
obtained with the corresponding single grid. However, as explained in the 
previous subchapter, there appears a narrow additional transmission win- 
dow due to interaction between the grids (this corresponds to the reflection 
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Figure 20: Reflection and transmission coefficient for a double grid of wires 
loaded with a parallel connection of bulk capacitors and inductors. 
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(a) Material parameters over a wide frequency (b) The same result as in (a), but plotted over a 
range. more narrow frequency range. 



Figure 21: Mesoscopic effective permittivity and permeability of a double 
grid of wires loaded with bulk a parallel connection of bulk capacitors and 
inductors. 



minimum around d/X ~ 0.34 point). 

The mesoscopic effective material parameters are shown in Fig. [2TJ All the 
parameters are purely real and grow as a function of frequency outside the 
resonance regions. Qualitative description presented in connection with the 
corresponding single grid still holds, however, as is clear from the discussion 
in the previous subchapter, the out-of-phase currents flowing in the grids 
leads to a resonant effective permeability. Again over a narrow frequency 
band just above the permeability resonance both material parameters are 
negative. 

Fig. [22] compares the effective material parameters obtained using the 
proposed method and the S'-parameter retrieval procedure for the case when 
the slab thickness equals the physical thickness. When inspecting the pa- 
rameters obtained through the S-parameter retrieval procedure we observe 
the branch point of effective refractive index at d/X « 0.376. However, even 
when excluding the disturbance caused by the compensation of the branch 
jump the permeability function, obtained using the S'-parameter retrieval 
procedure, is clearly nonphysical. The reflection/transmission results for the 
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homogeneous slabs (the effective parameters are those shown in Fig. |2"TI) are 
shown in Fig. [23J The behavior is qualitatively similar to the behavior with 
a double grid of capacitively loaded wires (please see the discussion in the 
previous subchapter). 



4 Conclusion 



We have studied single and double grids of wires loaded with a certain dis- 
tributed reactive impedance as an example of thin composite slabs. The 
currents induced by a normally incident plane wave to all the wires have 
been rigorously calculated taking into account the full interaction between 
the wires. Using the aforementioned data we have calculated the average 
polarization and magnetization induced to the structures, and the averaged 
fields confined in the unit volume. Using these quantities we have intro- 
duced mesoscopic effective material parameters for the proposed thin com- 
posite slabs. A large set of calculated example results has been provided. 
The physical background of the introduced mesoscopic material parameters 
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(a) Comparing the mesoscopic effective permittiv- (b) Comparing the mesoscopic effective permeabil- 
ities, ities. 

Figure 22: Comparison between the mesoscopic effective material parameters 
obtained using the proposed method and the S -parameter retrieval method. 
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Figure 23: Comparison between the reflection data calculated for the actual 
double grid of LC loaded wires and for the corresponding homogenized slabs. 



has been discussed in details, and the proposed results for the double grid 
structures have been compared with the results predicted by the commonly 
adopted ^-parameter retrieval procedure. 
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